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Abstract 



We analyze the g 2 -dependence of the strange magnetic form factor, 
Gm (g 2 ), using heavy baryon chiral perturbation theory (HBxPT) and dis- 
persion relations. We find that in HBxPT a significant cancellation occurs 
between the 0{p 2 ) and 0(p 3 ) loop contributions. Consequently the slope 
of G^m a t the origin displays an enhanced sensitivity to an unknown 0(p 3 ) 
low-energy constant. Using dispersion theory we estimate the magnitude of 
this constant, show that it may have a natural size, and conclude that the 
low-q 2 behavior of G$ could be dominated by nonperturbative physics. We 
also discuss the implications for the interpretation of parity- violating electron 
scattering measurements used to measure (g 2 ). 
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The flavor structure of low-energy nucleon properties (e.g., mass, spin, etc.) continues 
to be a subject of considerable interest in hadron physics. In particular, evidence exists 
that strange quarks may play a non-negligible role in such low-energy properties [[!]]. The 
strange quark contribution to the nucleon's electromagnetic structure is under investigation 
using parity- violating electron scattering experiments at MIT-Bates 0, the Jefferson Lab 
(JLab) [@,f|, and Mainz |J. These experiments have produced the first determinations of 
the strange magnetic form factor, Gm (g 2 ) 0, as well as a linear combination of Gm (g 2 ) 
and the strange electric form factor, Ge\q 2 ) JU- 

Theoretical interest has focused largely on the leading moments of these form factors: 
the strange magnetic moment, ji s = Gm\o) and the corresponding electric and magnetic 
radii, (v 2 )e,m = GdGE^ M (q 2 ) /dq 2 at q 2 = (note that we do not normalize the radii to the 
q 2 = values of the form factors). While the radii are intrinsically interesting, knowledge of 
(v 2 )m is also needed in order to extract fi s from experiment, since Gm can be determined 
experimentally only for non-zero q 2 . In order to extrapolate G$ to the photon point, one 
might hope to rely on chiral perturbation theory, which is in principle well-suited to charac- 
terizing the leading g 2 -dependence of form factors. While (v 2 )m is nominally of chiral order 
p 3 (we count the magnetic moment operator as being 0(p)) ) the leading order contribution 
from kaon loops is proportional to l/mx ~ and, thus, of 0(p 2 ). At this order, 

there exist no analytic (in quark mass) operators and, consequently, no counterterm or cor- 
responding low-energy constant. Thus, as observed by the authors of Ref. |J, one is able to 
make a parameter- free prediction for (r 2 )M at this order which may be used to perform a 
model-independent extrapolation of G$ to the photon point. This procedure has been used 
to extract a value of /i s = 0.01 ±0.29 ±0.31 ±0.07 from the value of the form factor obtained 
by the SAMPLE Collaboration at q 2 = -0.1 (GeV/c) 2 : G { $(q 2 = -0.1) = 0.14±0.29±0.31 
(The last error in fi s corresponds to the theoretical extrapolation uncertainty quoted in 
Ref. §.) 

Here, we compute the 0(p 3 ) loop contributions to G^iq 2 ) in heavy baryon chiral per- 
turbation theory (HBxPT). We show that the 0(p 3 ) contributions to {t 2 s )m largely cancel 
the 0(p 2 ) term. In addition, we observe that a new magnetic radius operator arises at this 
order, whose coefficient cannot presently be determined apart from a determination of G$ 
itself. Because of the cancellation between the 0(p 2 ) and 0(p 3 ) terms, the relative impor- 
tance of this coefficient, or low-energy constant (LEC), is enhanced over what one might 
otherwise expect. Using a dispersion theory analysis of Gm , we estimate this LEC and argue 
that it could give the dominant contribution to {t 2 s )m- As a corollary, we suggest that the 
low-momentum behavior of G$ is governed by non-perturbative physics and observe that 
one cannot presently extract fi s from experiment in a model-independent manner. 

In HB^PT, the contributions to Gm (q 2 ) through 0(p 3 ) are generated by the graphs 
of Fig. 1. The 0(p) contributions to fi s arise from insertion of the lowest-order magnetic 
operator in Fig. 1(a). The 0(p 2 ) contributions nominally arise from Figs, l(b-e). However, 
graphs l(c-e) contribute only to the charge form factor, so that the 0(p 2 ) contributions to 
fi s arise entirely from Fig. 1(b). The 0(p 3 ) contributions to fi s arise from several sources: 
insertion of the tree-level magnetic moment operator in Figs. l(f,g); operators containing 
two-derivatives (Fig. 1(h)); 1/M/y, or "recoil", corrections in Figs, l(i-k) (denoted by the 
a x"); and magnetic operators proportional to m s (Fig. 1(1)). We note that the sum of 
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Figs. l(j) and l(k) vanishes for Gm . 

Similarly, the leading contribution to the slope of Gm (g 2 ) at the origin also arises from 
Fig. 1(b) and is 0{p 2 ). The 0(p 3 ) contributions to the slope are given by the sum of Fig. 1(h) 
and l(i), along with the insertion of the magnetic radius operator in Fig. 1(a). One must 
also include the dependence of Fig. 1(b) on v • q = q 2 /Mat- 

The formalism for evaluating these diagrams in HB^PT has been discussed extensively 
elsewhere, so we do not give all the relevant formulae here (see, e.g., Ref. [|7| and references 
therein). However, in order to define our normalization for the LEC's we give expressions 
for the relevant magnetic moment and radius operators: 

A = ^VpXTr (BS°B) Z*" (1) 



W ^{6+Tr (BS a {Q, B}) + 6_Tr (SS a [Q, B 



A x 



Csb = W ^F^{& 3 Tr (bS°[[Q, B],M 



+b 4 Tr(BS a {[Q,B],M})+b 5 Tr(BS a [{Q,B},M\) (2) 
+b 6 Tv (BS°{{Q, B}, M}) + b 7 Tr (bSTB) Tr (MQ)} 

+^ W X^Tr (BS°B) Tr (M) 
£2 = -rJW</Tr (bS°b) d 2 Z^ (3) 

X 

W ^{& r + Tr (BS°{Q, B}) + b r _Tr (SB*®, B]) }d 2 F^ 

A x 

where B denotes the baryon field octet, Q = diag(2/3, — 1/3, — 1/3) (at lowest order), 
Ai = diag(0, 0, 1) is the SU(3) symmetry-breaking matrix associated with the light quark 
masses, v p and S a are the baryon velocity and spin, respectively, F^ v is the photon field 
strength tensor, and Z pv is the corresponding tensor for a source coupling to baryon number 
current 1 . The chiral scale A x = 4^^, with F n = 92 A MeV. The full chiral structure of the 
electromagnetic charge operator is obtained by making the replacement Q —> ^(^Q£+£Q^) : 
where £ = exp(in/F 7r ) and fl is the octet of pseudscalar mesons. Note that the symmetry- 
breaking operators in Csb are suppressed by two powers of p relative to C\. Following 
standard conventions, however, we absorb this suppression into the definition of the con- 
stants 63-8- Indeed, when b± and the octet constants bs-7 are determined from fits to octet 
baryon magnetic moments [0] , the latter are generally suppressed by an order of magnitude 
relative to the former, as one would expect from the p 2 suppression of Csb- 
The two-derivative operators have the structure 

Cmb = (BS°A») Tr (A"B) (4) 

A X 



1 A typo appears in Eq. (17) of Ref. Q. The term proportional to b& should contain a double 
ant i- commut at or . 
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-6i Tr (BS°[A», A v }B) + b u Tr [BS a {A^ A v }B)} , 



where A^ = (i/2)(£t<9^ - fd^*). 

Using these normalizations, and expanding 



GM(g 2 )=/, s + ^<r s W + "-- , (5) 



we obtain 




2 2 

m v , mi 



4.(1+ E (7 m -AW)^ln^f] (6) 



+ (5D 2 - 6DF + 9F 2 ) f ^ - In 41 

+b 8 - 2(b 3 + b 4 - \b, - \b, - b 7 ) - h 7 ^ In 

<^ 2 >m = ~< ( ^ ) K + ^ - 6DF + 9F>) g + Tin ^) (7) 
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where 

6 S = b - 2[6_ - (6 + /3)] 

jr = &S _ 2 [fir - (6;/3)] (8) 
h = b 9 — 26io + 66n 

and where /i is the subtraction scale associated with the divergence of the Feynman diagrams, 
D and F are the SU(3) axial current reduced matrix elements, and 7^ and \( x > are known 
functions of D and F JJJ. The dependence of b s (Jf s ) on the SU(3) singlet constant b (6q) 
and octet constants 6± (&±) 3 as we U as the presence of the 6 3 _ 8 terms in fi s arise from the 
group structure of the strange quark vector current: 

s 1 ,s = J*-2J I =\EM) , (9) 

where is the SU(3) singlet, or baryon number, current and J^ ={) (EM) is the isoscalar 
electromagnetic (EM) current, an SU(3) octet operator. Nucleon matrix elements of the 
latter receive contributions from the 6 3 _7 operators in Eq. (^), as evaluated in Ref. [^]. The 
operators proportional to 6 ,8 and % contribute to (p\J^\p). While values for the octet LEC's 
&± 3 b±, and 6 3 _7 can be obtained from measured octet baryon magnetic moments, one has 
no independent determination of the singlet LEC's 6 , 6q, and 6 8 (apart from measurements 
of Gm itself). Note also that the LEC's are functions of the subtraction scale, bi = bi{fi) 
and b s = 6 5 (/x). 

The contributions through 0(p 2 ) for both ji s and (r^M have been computed in Refs. |||8|. 
In the case of (r^, the 0{p 2 ) contribution contains the factor M N /rriK> The contributions 
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of 0(p 3 ) have not been computed previously. As a check on our results, we compare with the 
integrals appearing in the relativistic calculation of the isovector electromagnetic form factors 
in Ref. |J. After replacing by m^, we find agreement with the non- analytic contributions 
through 0(p 3 ). Note that symmetry-breaking effects which break the degeneracy of the 
pseudoscalar decay constants (F n ^ F K ^ F v ) contribute at higher order than we consider 
here, so we use only in our numerical evaluation below. Similarly, the lowest-order effect 
of the Ma,e and M N non-degeneracy is analytic in m s and can be absorbed entirely in the 
LEC's &o, b±, &o> e ^ c " Finally, as in Ref. [|], we have included only octet baryon intermediate 
states. While inclusion of decuplet states will likely affect the precise numerical values of the 
0{p 2 ) and 0(p 3 ) contributions, we do not expect modification of our qualitative conclusions. 

The low-|g 2 | behavior of G$ depends on several LEC's, including 6 S , b r s , and frs-n- To 
determine the relative importance of the LEC's and loops, we take D = 0.75 and F = 0.5, 
set ji = 1 GeV, and derive values for 6 3 . . . 6 7 from the octet magnetic moments. In doing 



so, we evaluate the combination of constants h as in Ref. |10| . The result gives 



\l s = 2.5 [b a (fi = 1 GeV) + 0.66 8 (/x = 1 GeV)] + 1.2 (10) 
{^ s )m = ~ [0.04 + 0.3^(/i = 1 GeV)] fm 2 , (11) 

where the LEC-independent contributions 1.2 (in /x s ) and —0.04 fm 2 (in {t 2 s )m) arise entirely 
from the loop graphs and symmetry-breaking LEC's. We note in particular the relative 
sensitivity of (t 2 s )m to loops and the LEC b r s . Were the latter to have a "natural" size, 
\b r s \ rsj l ? its contribution would dominate the slope of Gm and the origin. Indeed, an 
uncertainty in this LEC of Ab r s = ±1 corresponds to a ±0.3 fm 2 uncertainty in the radius. 
This situation contrasts starkly with the result at 0{p 2 ) ) which gives (t 2 s )m = —0.16 fm 2 
(arising entirely from Fig. 1(b)). In short, inclusion of the 0(p 3 ) loop contribution produces 
a substantial cancellation with the 0(p 2 ) term, thereby enhancing the relative sensitivity of 
( r s)M to the LEC b r s . While one generally expects a suppression of higher-order effects, the 
rate of convergence for the strangeness form factors is governed by m^/A x ~ 1/2 and, thus, 
is rather slow. Consequently, the occurance of significant cancellations between successive 
orders as we find here is not entirely surprising. 

At this time, it is not possible to determine ji s and (t 2 )m from existing measurements. 
As evident from Eq. (|9|) the strange quark vector current contains both SU(3)-octet and 
SU(3)-singlet components. Contributions to to b s + 0.66s and b r s from the octet component 
are known, while the singlet component is undetermined. In principle, one might attempt to 
determine /i s - and, thus, b s + 0.66 8 - from the published value of G$(q 2 = —0.1), using the 
lowest-order HB^PT result for (v 2 )m as was done in Refs. However, as the results in 

Eqs. (|7[,|TTD indicate, such a procedure could be misleading. Given the sensitivity of (t 2 )m to 
b r s) neither the magnitude nor the sign of the magnetic radius can be determined at present 
in a model-independent manner. Thus, the correct extrapolation of Gm from q 2 = —0.1 
(GeV/c) 2 to the photon point is simply not known at present. 

To illustrate this situation, we plot in Fig. 2 G$ as a function of Q 2 = — g 2 , showing 
the published result at Q 2 = 0.1 (GeV/c) 2 . We also show various possible extrapolations 
to Q 2 = 0. The dashed lines indicate the extrapolation for b r s = — 1 and the dot-dashed 
lines give the extrapolation for b r s = +1. Taking these two possibilities as reasonable, 
though not rigorous, extremes, one would infer that the strange magnetic moment falls in 
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the range 0.7 > fi s > —0.4. Using only the leading order HB^PT results one would infer 
0.5 > /i s > —0.3. We emphasize, however, that larger values for \b r s \ than used in Fig. 2 
would not be unreasonable, so we would not want to use Ab r s = ±1 as an estimate of the 
theoretical uncertainty. For example, changing b r s from ±1 to ±2 increases the range in ji s 
by ~ ±0.1. We also note that, while the present uncertainty in ji s is likely dominated by 
the experimental error, the relative importance of the |&g|-induced uncertainty could become 
comparable to the experimental error with more precise data for G^iq 2 = -0.1). 

In general, LEC's such as b s and b r s characterize the effects of short-distance physics 
arising from scales r < 1/A X . Such non-perturbative effects can, in principle, be computed 
using lattice QCD or dispersion relations. While no definitive calculation has been performed 
using either approach, it is nevertheless instructive to compare the implications for the 
LEC's from existing analyses. In the quenched lattice calculation of Ref. JTT[] , Gm (g 2 ) was 



evaluated at five different kinematic points, starting at Q 2 = and increasing by increments 
of roughly AQ 2 = 0.4 (GeV/c) 2 . The results at Q 2 = give ji s = — 0.36±0.20, corresponding 
to b s + 0.6&8 ~ —0.6 ±0.1. Given both the spacing of the lattice data as well as the size 
of the error bars, there exists considerable lattitude for the shape of G$ at low-Q 2 . For 



example, the authors of Ref. [|TT|] fit the lattice data to a dipole form factor, yielding a radius 
parameter = 1 GeV) = 0.23 ± 0.13. However, a different shape for G$ a t low-Q 2 - 
including an opposite sign slope at the origin - would not be inconsistent with the published 
lattice results. 

In the case of dispersion relations, the strange magnetic moment and magnetic radius 
are related to integrals over the imaginary parts of the relevant strange form factors: 

M* = - / 12 

r s ) M - - I dt . (13) 
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Note that since the Sachs form factor Gm \q 2 ) is the sum of Dirac and Pauli form factors, 
F[ s ^ and F 2 , respectively, and since f[ s \q) = 0, only the integral over contributes to 

Ms- 

In each case, the integrands receive contributions from a tower of hadronic states: 37T, 
57T, 77T, KK ) .... In previous work |1^,[T3|], we evaluated the KK contribution using KN 



scattering and e + e~ —> KK data. We find that the a kaon cloud" contribution to both the 
strange vector form factors and EM isoscalar form factors is dominated by the 0(1020)- 
resonance. For the purely pionic intermediate states, e + e~ data suggest that only the 3tt 
channel gives an important contribution and that it is dominated by the cj-resonance. Using 



the known flavor content of the a;, one can follow Ref. [|T4| and evaluate its contribution to 

Ms) 

• 

Inclusion of higher-mass states (e.g., KKtv) is necessary in order to obtain a physically 
reasonable large-Q 2 behavior of vector form factors. In the case of the isoscalar magnetic 
form factor, the sum of the 3n and KK contributions - along with a single higher mass 
pole to approximate the effect of the remaining higher-mass states - gives an excellent 
description over a wide range for Q 2 . Since we do not compute the higher-mass contributions 
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from first principles, however, we are not confident in evaluating their contribution to . 
Consequently, we quote only the "low-mass" (3tt and KK) contributions to the moments 
in Eqs. (0, [L3|): fi s = -0.28 and (r 2 s ) M = -0.34 fm 2 0. These values correspond to 
b s + 0.6&8 ~ —0.6 and 6^ « — 1.1, respectively. A comparison of these values with the 
lattice/dipole fit values is given in Table I. The corresponding low-mass dispersion relation 
prediction for the low-Q 2 behavior of G$ is also shown in Fig. 2. A reasonable, though 
not rigorous, estimate for the uncertainty associated with presently unknown higher-mass 
contributions is also shown (shaded band)0. For example, a negative slope but with smaller 



magnitude arises in the pole approximation analyses of form factors given in Refs. [Tj,[T5[ . 

We observe that the lattice and "low-mass" dispersion relation determinations of b s are 
comparable, while the dipole fit to lattice data gives a significantly different prediction for the 
radius parameter b r s than obtained with the dispersion relation analysis. We speculate that 
this difference could be an artifact choosing a dipole form factor to parameterize the lattice 
evaluation of Gm\ a result of the omission of higher-mass contributions to the dispersion 
integral, or both. More importantly, each of these approaches, which incorporate both chiral 
and nonperturbative effects, produce values for the LEC b r s of natural size, though having 
opposite sign. In either case, one would conclude that non-perturbative dynamics play an 
important role in the physics of Gm • Ideally, future theoretical progress will ultimately 
lead to consistency between the lattice and dispersion relation treatments of these non- 
perturbative effects. Such progress will undoubtedly require incorporation of the correct 
chiral structure of G$ at low-|g 2 |. From the experimental side, additional data for Gm 
at larger \q 2 \ may not be sufficient to determine fi s and (r 2 )M, since the chiral expansion 
is valid for a limited |g 2 |-range (e.g, out to ~ m 2 K ) and since new operators contribute as 
one moves successively further from the photon point. On the other hand, additional low- 
\q 2 \ data [e.g., \q 2 \ < 0.05 (GeV/c) 2 ] from parity-violating electron scattering experiments 
could provide an unambiguous determination of the behavior of strange magnetism at small 
momentum-transfer. 

This work was supported in part under U.S. Department of Energy contract #DE-FG02- 
00ER41146 and by the National Science Foundation under award PHY00-71856 and grant 
No. PHY-0098645. We thank B.R. Holstein, B. Kubis, R.D. McKeown, U.-G. MeiBner, and 
T.M. Ito for useful discussions. 



2 Additional potential sources of uncertainty include those associated with the cj-pole approxi- 
mation for the 3tv contribution, the analytic continuation used to obtain the spectral functions, 
etc.. 
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FIG. 1. Contributions to strange magnetic moment and radius. Circle denotes leading-order 
magnetic moment and radius operators, while square indicates operators proportional to quark 
mass or containing two derivatives. The x indicates a 1/Mjy insertion. Solid line and dashed lines 
indicate octet baryons and pseudoscalar mesons, respectively. 




Q 2 [GeV 2 ] 

FIG. 2. The Q 2 -dependence of Gm . Circle at Q 2 = 0.1 (GeV/c) 2 indicates the result from the 
SAMPLE measurement [Q. Inner error bar gives statistical error while outer error bar combines 
statistical and systematic errors in quadrature. Dashed (dot-dashed) lines indicate extrapolation to 
Q 2 — using radius parameter b r s — —1 (+1). Solid line gives "low-mass" dispersion relation result 
fl3[ , while shaded region indicates possible effects of higher mass contributions to the dispersion 
relation. 
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TABLES 



Method 


b s ± 0.6& 8 


K 


Lattice/Dipole Fit 


-0.6 ±0.1 


0.23 ±0.13 


Low-mass Dispersion Relation 


-0.6 


-1.1 



TABLE I. Low-energy constants for strange magnetic moment (b s ) and strange magnetic 
radius (b r s ) obtained from (a) lattice QCD calculation of Ref. fll| and (b) 3tt and KK contributions 
using dispersion relations [Jl3| . The the lattice value for b r s was obtained by fitting lattice data to 
a dipole form for Gm (g 2 ). 
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